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Abstract 

We compute the gap probability that a circle of radius r around the origin contains exactly 
k complex eigenvalues. Four different ensembles of random matrices are considered: the Ginibrc 
ensembles and their chiral complex counterparts, with both complex (/3 = 2) or quaternion real 
(f3 = 4) matrix elements. For general non-Gaussian weights we give a Fredholm determinant or 
Pfafhan representation respectively, depending on the non-Hermiticity parameter. At maximal non- 
Hermiticity, that is for rotationally invariant weights, the product of Fredholm eigenvalues for f3 = 4 
follows from /3 = 2 by skipping every second factor, in contrast to the known relation for Hermitian 
ensembles. On additionally choosing Gaussian weights we give new explicit expressions for the 
Fredholm eigenvalues in the chiral case, in terms of Bcsscl- K and incomplete Bessel-I functions. 
This compares to known results for the Ginibre ensembles in terms of incomplete exponentials. 
Furthermore we present an asymptotic expansion of the logarithm of the gap probability for large 
argument r at large N in all four ensembles, up to including the third order linear term. We can 
provide strict upper and lower bounds and present numerical evidence for its conjectured values, 
depending on the number of exact zero eigenvalues in the chiral ensembles. For the Ginibrc ensemble 
at (3 = 2 exact results were previously derived by Forrester. 



1 Introduction 



Non-Hermitian Random Matrix Theory (RMT) as introduced by Ginibre pQ is almost as old as the 
classical Wigner-Dyson ensembles. In the past decade it has enjoyed a revival of interest and many 
different applications of it have been made, where we refer to [2] for a recent review. However, most 
works have concentrated on computing the spectral correlation functions, culminating recently in the 
solution of the Ginibre ensembles of real asymmetric matrices [3] . 

In this article we will focus on the computation of gap probabilities from which the distribution 
of individual eigenvalues or spacings can be derived. Once the spectral correlations of an ensemble 
are known in terms of a kernel of (skew) orthogonal polynomials in the complex plane, one can in 
principle express the gap probabilities in terms of a Fredholm determinant involving this kernel, as we 
will show. However, such expressions are in general not very explicit, unless the Fredholm eigenvalues 
are known. On the other hand an expansion of the Fredholm determinant in terms of integrals 
over spectral correlation functions converges very rapidly (see e.g. [1]), but nevertheless becomes 
cumbersome because of having 2D integrals in the complex plane. 

Explicit results for the gap probability were first derived for the (3 = 2 Ginibre ensemble with 
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unitary invariance at maximal non-Hermiticity [5]. Here the probability Eq (r) that a circle of radius 
r around the origin is empty is given in terms of a product over incomplete exponentials. The same 
quantity follows for the (3 = 4 Ginibre class [6], where every second term in the product is skipped. 
This result was then used in [5] for (3 = 2 to compute the so-called level spacing distribution in the 
complex plane, by placing one eigenvalue at the origin and computing the probability to find a second 
eigenvalue at radius r. Due to the translational invariance with respect to the large-iV macroscopic 
density being constant on a disc the corresponding repulsion of complex eigenvalues is supposed to 
hold everywhere in the bulk [5]. 

The first aim of the present work is to extend the above results to the chiral complex ensembles at 
(3 = 2 and 4, computing explicitly the product of Fredholm eigenvalues at finite-iV. These ensembles 
were introduced in [7] {(3 = 2) and [8] {(3 = 4). One motivation for this is the application of non- 
Hermitian RMT to Quantum Chromodynamics (QCD) with non-zero chemical potential, and we refer 
to [9] for references and a review. The RMT predictions in the complex plane can be compared to 
numerical solutions of QCD. The first comparison on C was done using the complex level spacing 
distribution of [5] of the (3 = 2 Ginibre ensemble in the bulk of the spectrum [10]. However, in QCD 
chiral symmetry is very important at the origin, and a successful comparison to individual complex 
eigenvalues there was made very recently [3]. Our results for the gap in the chiral complex (3 = 2 
class valid at maximal non-Hermiticity were already announced there. We will give a derivation of 
this result and extend it to (3 = 4. 

Apart from these explicit results we give a Fredholm determinant (Pfaffian) form valid for a general 
weight function and non-Hermiticity parameter at (3 = 2 (4). This is shown to be equivalent to a 
(matrix) eigenvalue equation involving the kernel of (skew) orthogonal polynomials. For rotationally 
invariant weights it implies that the relation between the gap at (3 = 2 and 4 skipping every second 
Fredholm eigenvalue holds in general, extending this relation known for the Gaussian Ginibre case [6]. 

The second part is then devoted to the asymptotic expansion of the logarithm of Eq{t) in all 
four Gaussian ensembles at maximal non-Hermiticity. Our motivation here is to extend similar con- 
siderations done for Hermitian RMT. These include the so-called Widom-Dyson constant to the third 
order, and the first computation in [11] was made rigorous only very recently |12j . 

Our approach will be more heuristic in the sense that we give an exact derivation of the first two 
terms using standard asymptotic techniques as the Stirling formula. For our corresponding third order 
linear terms we can provide strict upper and lower bounds for their coefficients, as well as numerical 
evidence for their conjectured values. Our results are consistent with the exact results for the Ginibre 
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ensemble at j3 = 2 given in [13] to the fourth order. The chiral ensembles which are our main focus 
depend explicitly on the number of exact zero-eigenvalues in the third order linear term. 

Our paper is organised as follows. In Section [2] we define our ensembles and correlation functions. 
Section [3] first introduces Fredholm theory in the generic non-rotationally invariant case, followed by 
a calculation of the new product formulae for E^\r) for both Gaussian chiral complex ensembles at 
maximal non-Hermiticity. Section H] is devoted to a detailed analysis of our asymptotic expansion of 
log[E^f\r)] for Ginibre /? = 2 at large N, and Section [5] summarises and compares the corresponding 
results for the other three ensembles, including a small radius expansion. Our concluding remarks in 
Section [6] are followed by several appendices where some technical results are collected. 



2 Definitions 

In this section we collect together the definitions of the random matrix ensembles to be considered, 
starting with the Gaussian Ginibre ensembles and then moving to their chiral complex counterparts. 
We then define all density correlation functions, gap probabilities and individual eigenvalue distribu- 
tions and point out their mutual relationships. 



2.1 Ginibre ensembles 



The Ginibre ensembles depending on a non-Hermiticity parameter r, also known as Ginibre-Girko or 
elliptic ensembles, are defined as (see e.g. [HI [IS]) 



Z Gin ~ 



d& exp 



1-r 2 V 2 V 



, re [0,1) . 



(2.1) 



Here $ is a complex non-Hermitian N x N matrix for /3 = 2 or quaternion real matrix for (3 = 4. This 
ensemble can also be thought of as a two-matrix model $ = H\ + WH2, with v = [(1 — t)/(1 + r)]5 , 
being composed of two Hermitian (self dual) matrices H\ t i with distribution exp[— NTrHj/(l + r)]. 
In the limit r — ► l the Gaussian Unitary or Symplectic Ensemble is recovered. 

In both ensembles one can go to a complex eigenvalues basis Zj=i jv of the matrix $ 



2, 



08) 
Gin 



N 

n / d 2 zj WGin ( Zj ) jg\{z}) 

j=l J C 



(2.2) 



with Jacobian 

r(/?=4) 



H\z k - Zl \ 2 = \A N ({z})\ 2 , 

k>i 

N N 



N 



J, 



N d Z }) = Yl\ z k~ Zl\ 2 \ Zk -Zi\ 2 
k>l 



n \ zh - z h\ 2 



h=l 



Y[(z h - z* h )A 2N ({z,z*}) 
h=l 



(2.3) 
(2.4) 



Both can be expressed through a Vandermonde determinant Ajv(z), but in a way different from the 

Hermitian ensembles with real eigenvalues. In A2n({ z , z*}) the variables are ordered as z%, zf, z 2 , 

The exponential weight function is given simply by 



WGin{z) 



cxp 



N 
1^72 



^+z* 2 ) 



(2.5) 
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for both (3 = 2, 4. In the case r = called maximal non-Hermiticity it becomes rotationally invariant, 
a case we will study in great detail later. 

At P = 2 the ensemble can be solved using Hermite polynomials as orthogonal polynomials in the 
complex plane, and we refer to [16] for an exhaustive discussion. At j3 = 4 the solution is given in terms 
of skew orthogonal polynomials in the complex plane, again constructed from Hermite polynomials, 
and a complete discussion is given in [T7]. At r = a basis of (skew) orthogonal polynomials is always 
constructed from monic powers instead of Hermite polynomials. 



2.2 Chiral ensembles 

Next we turn to the chiral complex counterparts of these Ginibre ensembles. They are defined as 
Gaussian two-matrix models [8] 



d^d^ exp 



4/i 2 



Tr(^$ + $t^t) 



fi E (0, 1] . 
(2.6) 

Here $ and ^ are two matrices of rectangular size (N + v) x N with either complex (/? = 2) or 
quaternion real {(3 = 4) elements, without further symmetry properties. The two matrices can be 
composed from two matrices $ = iH\ + [iB.2-, and ^ = iH\ + f-iHr, with Hi2 being non-Hermitian 
(quaternion real) with distribution exp[— NTrHjHj] each. In the limit /i — > leading to $ = — \pt we 
recover the chiral Gaussian Unitary or Symplectic Ensemble 0. 

We are interested in the complex eigenvalues of the matrix ^ ^ ^, or equivalently of 

This change of variables is tedious, by first going to complex eigenvalues of $ and each and then 
integrating out one set. With details given in [7j ((3 = 2) and [8] {(3 = 4) we only give the answer 



with the following non-Gaussian weight function 



N 
7 i 



(2.7) 
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K 



N\z\ 



2/i 2 



cxp 



N(z z + z* 2 ) 



1 



4// 2 



(2.8 



which now depends on (3 and the number of exact zero-eigenvalues v > 0. We note that the Jacobian 
i7 only differs from the Ginibre ensembles by inserting squared variables. 

The chiral ensembles eq. (|2.7p can be solved in terms of Laguerre polynomials as orthogonal 
[7j [T8l E] and skew orthogonal [8] polynomials in the complex plane for (3 = 2 and 4 respectively. 

Regarding the weight we note that the Bessel-.fr function again becomes an exponential at half 



integer values of the index, e.g. K, i(a 



lirn^oo K u (x), or in the asymptotic limit of 



large arguments which is also reached when taking \i — > 0. In the opposite limit \i 
non-Hermiticity we obtain again a rotationally invariant weight. 



1 at maximal 



2.3 Correlation functions 

Let us now turn to the correlation functions to be calculated, which we define for arbitrary weight 
functions. We define the fc-th gap probability with respect to radial ordering to give the probability 

lr The non-Hermiticity parameters of the two sets of ensembles can be brought onto an equal footing by mapping 
M 2 = (1-t)/(1 + t). 
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that k (independent) complex eigenvalues lie inside the circle of radius r around the origin, and N — k 
lie outside. Here we do not count those v eigenvalues which are always zero by construction in the 
chiral ensembles. For the chiral ensembles the definition reads El, 

1 Ap k „ r N N „2-7r 

Ef(r) = Tm (N _ kW U dr.r, J] / dr in \{ d0 n wW(z n ) jtf\{z*}) , (2.9) 

Z ch { ' j=l J ° l=k+l Jr n=l J ° 

where we have switched to polar coordinates z n = r n e tdn . The quantities for the Ginibre ensembles are 
the same and simply obtained by replacing the last two factors by WGin{z n )Jiq\{z}), the respective 
joint probability distribution function (jpdf). 

In eq. (|2.9|) we could have given a more general definition, by choosing other sets (i.e. not 
necessarily concentric circles around the origin) or by including an angular dependence. While this 
was briefly sketched in [5] we focus here on quantities that we can compute most explicitly. We will 
come back to this point in the next section. 

From the fc-th gap probability defined above, the radial distribution of the k-th individual eigen- 
value ordered with respect to radius follows recursively by a simple differentiation 

JUf (r) = r (pf(r) - pjg^r)) , (2.10) 

ix this ir 

relation eq. (|2.10p is easily inverted to 



setting Pq\t) = 0. In particular this implies for the first eigenvalue that -t^Eq (r) = —rpf\r). The 



1 d k 



r dr 

1=0 



Obviously it should hold that the sum over all individual eigenvalues gives the spectral density (z) 
to be defined below, 

/*27T f OO 

/ d6Rf ] (z = re ie ) = Vpf } (r) , with / drrpf\r) = 1 Vfc . (2.12) 
Jo V / k=1 Jo 

Here we have integrated over the angle as the Pu\r) only depend on the radius. They are normalised 
with the two-dimensional radial measure. This relation will be used later to illustrate and check our 
results for individual eigenvalues. We now come to all the Appoint density correlation functions defined 
as 

i m\ N r N 
z W(N-k)\ .^ +1 J C f = \ 

with the same quantities for Ginibre defined using their respective jpdf. In particular for k = 1 the 
spectral density is normalised to the number of eigenvalues: f c d 2 zR^f\z) = N. 

All fc-point density functions are known explicitly for finite iV for all four ensembles defined above: 

Rf\ Zl ,...,z k ) = (Q)det 1 < M .< fc [Kjf ) (z J ,z*)] , (2.14) 

(2) 

where K^'(zi, z*) is the kernel of the corresponding orthogonal polynomials in the complex plane for 
(3 = 2. In our Gaussian ensembles these are of Hermite |15] type for Ginibre, or Laguerre I18j 



2 This definition of E ( k P) (r) differs from pjj] by a factor l/k\. 
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for the chiral case. For (3 = 4 the quaternion determinant Qdet (or Pfaffian) is taken over the 2x2 
matrix kernel expressed in terms of the pre- kernel KN(zi, zj) of the corresponding skew orthogonal 
polynomials in the complex plane. For our examples these are given again by Hermite [T7] (Ginibre) 
or Laguerre [5] (chiral) polynomials. 

The k-th. gap probability can be expressed in terms of n-point density correlation functions as 
follows 

E k\ r ) = h T. { - = 7rW F dr ^ I Rk+t(zi,...,Zk+i) , (2.15) 
e=o l - j=i Jo Jo 

as was pointed out in [J] for the chiral complex ensembles at f3 = 2. The term in the sum with k+£ = 
is set to unity. In fact this expansion holds for all 4 ensembles both at (3 = 2 and 4, independent of 
the structure of the jpdf (or even for real (3 in case of Hermitian RMT). This can easily be obtained 
using the same expansion as for real eigenvalues in [19], inserting merely the definition (|2.9p . and 
without using the solution eq. (|2.14l) . Moreover one can define a generating function by generalising 
the zero-th gap probability to 

N / poo rr \ f2ix 



leading to 



= E II / dr i r i / Mj ■ ■ ■ > H) , (2-16) 

£=0 ' j=l J ° 

Ef\r)= { -=^^ k E^{r-0\^ for k = 0,l,...,N. (2.17) 

The importance of this relation will become evident when we give a product representation in terms 
of the eigenvalues of a certain determinant. 

Note that although eq. (|2.15p is explicit with all terms on the right hand side known, it contains 
multiple integrals of determinants of an increasing size. As an approximation the sum can be truncated 
after the first few terms, and it was found in [1] to converge rapidly for k = 0. However, even within 
this approximation the few integrals to be done numerically become rapidly cumbersome. 



3 Results for the gap probabilities 

In Subsections 13.11 and 13 . 21 we give some general results for E^(r) valid for generic non-Hermiticity. In 
the chiral case this includes a determinant or Pfaffian formula for recursively known one-dimensional 
integrals. Then in Subsections l3.3l and l3.4l we turn to maximal non-Hermiticity where the determinants 
or Pfaffians can be diagonalised and an explicit product representation is given for any matrix size. 



3.1 General case (3 = 2 



We will start the discussion with j3 = 2 for a general weight function in the eigenvalues representation, 

(2) 

including the chiral and Ginibre case. First the gap probability Eq (r) is treated, and then through 
the generating functional E^ 2 \r;^) all subsequent gaps. Our derivation is in complete analogy to the 
real case, see Section 6.3 [6], but we will repeat it here in a slightly modified version to prepare for 
(3 = 4 which is new. We have that 

N „ N 



ZWch ^ J C \Cr f = l JC\C r 



d 2 Zl 



det liflj 
l<j,k,<N 



(3-1) 
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where C r denotes the circle of radius r around the origin. While this will become important in 
Subsections 13.31 and 13.41 throughout this subsection we could choose the set to be general. The 
Ginibre case is trivially obtained by choosing the other weight and non-squared arguments inside the 
Jacobian. 

Here we have also used that the normalising partition function is given as follows in terms of the 
squared norms hj of the corresponding orthogonal polynomials pj (z), f d 2 z (z)pj(z)pk(z*) = hj5jk, 

N-l 

Z { 2 = N\ [J hj . (3.2) 

3=0 

(2) 1 ~ — 

Examples for the orthonormalised wave functions <fj(z) = w v {z)^h- 2 pj(z) are 

^ch(z) = w W( z )h]h](J^) , (3.3) 
<PjGin(z) = ^Gl(^hf 2 Hj(zM) , (3.4) 



for the chiral ensemble and Ginibre ensemble, respectively. Applying Gram's result, e.g. see Appendix 
A. 12. of [S], we obtain 



JV-l 

(2), 



E}r , {r)= det 

l<j,k,<N 



S jk - [ d 2 z pj-iiz^k^iz*) = rj (l-Af } ) . (3.5) 

JC r J j =0 



This determinant has to be diagonalised. Although all integrals inside the matrix elements can be 
determined recursively for the chiral Gaussian case, see Appendix [A] for a slightly different framework, 

(2) 

this does not allow us to compute the eigenvalues A;- explicitly. Alternatively these eigenvalues can 
be obtained from the following integral equation 

Xif)(u) = [ d 2 zK^\u ) z*)^(z) , (3.6) 

JCr 

N-l 

K%\u,z*) = (3.7) 

3=0 

This is seen as follows. Because the l -fj{u) form a basis the eigenfunction ip(u) can be expanded in 
terms of them. Here, we will use a slightly more formal argument than in [6], by using the projection 
property of the kernel, 

r d 2 zK^\u,z*)K^(z,v*) = K^\u,v*) ,(3 = 2,4 . (3.8) 

As indicated this property will generalise to (5 = 4 (for the definition of the (3 = 4 kernel see eq. 
(pL20j) ). Integrating Kfi J (v, u*) times eq. (13. 6ft with respect to u and applying eq. (13. 8|) we have 

N-l r 
AV^)c» = / d 2 zK%\v,z*)^{z)=\^{v) (3.9) 



1=0 



with a = d 2 zifi(u*)tp(u) , (3.10) 
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and thus the desired expansion of the eigenfunction (for A / 0). Plugging this expansion back into 
eq. ()3,6p we have 

N-l N-l ~ N-l 

A Ci<p.i(u) = ^2 fj( u ) / d2z Pjiz*) ^2 Ci< ^( z ) • ( 3 - u ) 

■ n ■ n J Cr z n 



i=0 



j=0 



i=0 



The projection onto the coefficients q can be achieved using the orthonormality of the (fi(u), inte- 
grating both sides with f c d 2 uipi(u*): 



N-l . 

T=t Jc 



d 2 z ipi(z*)ipi(z) 



(3.12) 



This equation has a solution if 



= det 

l<i,j<N 



\5ij 



d 2 z ipi_ 1 (z*)(p j -i(z) 



(3.13) 



The determinant of this Hermitian matrix is given by a polynomial of degree N in term of its real 



eigenvalues X\ 



(2) 



det 

l<i,j<N 



X5i 



d 2 z tpi- X {z*)<Pj-i{z) 



N-l 



i=0 



n a-a 



(2) 



(3.14) 



This determinant is called the Fredholm determinant, and its eigenvalues are the Fredholm eigenvalues 
which obviously depend on r here. Setting A = 1 we are back to eq. ()3.5[) as we wanted to prove. 

Moving to the generating functional we can repeat all the steps. Multiplying all equations by £ 
and redefining A —> A£ it trivially follows for the generating functional 



JV-l 



i=0 



,(2) 



(3.15) 



Applying definition eq. (|2.16|) the following representation of the fc-th gap probability holds: 



N-l 



no-Af) e 



vG9) 

Si 



v09) 



j=0 



1 - A^ } 



.71 



i - \<a> 



(3 = 2,4 



(3.16) 



3k 



where we have anticipated that the same result holds for (3 = 4 (as will be shown next). The same 
result is equally true for Hermitian RMT, see [6] where it was derived in a different way. In eq. (|3, 161) 
the sum is over all possible permutations of subsets of k out of N indices, e.g. for k = 1 it reads: 



E?\r) = n(l-Af)E 



N-l 



N-l A (/3) 



3=0 



1-A 



09) • 



(3.17) 



3.2 General case (3 = 4 

Let us turn to j3 = 4. To the best of our knowledge the relation between the Fredholm eigenvalues and 
an eigenvalue equation involving the kernel is new. Looking at the Jacobian eq. f)2.4f) and replacing 
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the single Vandermonde by a determinant of monic, skew orthogonal polynomials and their complex 
conjugates (see eq. (j3.22j) ) we can write for the gap probability 



r(4) II 

' c h 1 = 1 
1 N 



Ml 



1 = 1 C\C r 



d 2 z l {z 2 -z* l 2 ) W M(z l ) 



det 

l<k<N, l<j<2N, 



(3.18) 



The 2N rows alternate in the variables Zk and zt- The iV-fold integral can be reduced to a Pfaffian 
over single integrals using the de Bruijn integration formula, see e.g in [17] . and we obtain 



^0 \ r > — —rli <k j <2 N 



z 



(i) 

ch 



d 2 z(z 2 - z* 2 ) W M(z)( Kqk . 1 (z)q l . 1 (z*) - q k ^(z*) qi ^(z)) 



. (3.19) 



This is the (3 = 4 result corresponding to eq. (|3.5p . valid for any weight function which we display 



explicitly here. The normalising partition function given by the squared norms, zj^ = (2N)\ IljlJo 1 
as can be seen from the same calculation but over the full complex plane, using the definition (|3.22p . 
The Ginibre ensembles are trivially obtained by choosing the corresponding weights and non-squared 
arguments inside the Jacobian and in front of the weight. 

Our next task is to relate the eigenvalues of the matrix inside the Pfaffian to eigenvalues involving 
the (3 = 4 matrix kernel given by [17j 



rJV-l 



K%\z,u*) ee (z*' z - z z )^{u*' z - u*)*wW{z)*w£>{u) 



* 2 



* 2 



,(4) ( 



,(4), 



where the pre-kernel is defined as 



kn(z*,u) —kn(z*,u*) 
kn(z,u) —kn(z,u*) 



N-l 



k n (z,u) = ^ h k 1 l92k+l(^)Q2k(u) - Q2k+i{u)q 2 k{z)] ■ 



(3.20) 



(3.21) 



fc=0 



The polynomials we choose to be skew orthogonal with respect to the following antisymmetric product: 



d z z(z* 2 - z z )w^ ] (z) q 2 k+i{z)q2l{z*) - q2k+i(z*)q2i(z) 



(3.22) 



and zero whenever two even or two odd polynomials are contracted. The kernel eq. (|3.20p defined in 
this way then satisfies the contraction property eq. (|3.8p [17] . 

We can now state the matrix eigenvalue equation that is needed to solve eq. (13.191) . 



A 



d z z(z* 2 - z z )w^{z) 



Cr 



(4), 



kn(u*,z) —kn(u*,z*) 
kn{u,z) —kn(u,z*) 



Hz*) 
Hz) 



(3.23) 



The two components are trivially related by complex conjugation and thus give rise to the same A. 
However, for the sequel it is useful to use this matrix form. As for (3 = 2 the generalisation to gaps 
with respect to other sets is trivial. 

To proceed the expansion of the eigenfunctions in analogy to eq. (|3.9[) can be achieved by mul- 
tiplying eq. (|3.23p by the matrix in eq. (|3.20p from the left, integrating over (u* 2 — u 2 )w^\u) and 
using the contraction eq. (|3.8p to arrive at 



A / d z u{u* 2 -u z )w^{u) 



(4), 



kn(v* ,u)ip(u*) — kn(v*, u*)ip(u) 
kn(v,u)tP(u*) — kn(v,u*)iP(u) 



A 



i(j(v) 



(3.24) 
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The desired expansion of the eigenfunction ip(v) in terms of the skew orthogonal polynomials thus 
reads 

N-i 1 

Tp(, v ) = ^ T-[l2k+i(v)c 2 k - q2k{v)c 2 k+i] (3.25) 



n h k 
k=0 



with 

Cfe 



d 2 u{u* 2 - u 2 )w^(u)[q k {u)^{u*) - q k (u*)4>(u)} . (3.26) 



Since the even skew orthogonal polynomials q2k( v ) are n °t uniquely determined (see [I?]), then neither 
is this expansion. However, after fixing the former there is no ambiguity in the definition of the c k s. 

Next we can project onto the coefficients q. Multiplying eq. (|3.23p with (—qi{u),qi(u*)) and 
integrating over u with the corresponding measure the skew orthogonality eq. (|3.22[) leads to 

„ N-l 
Aq = / d 2 z(z* 2 - z 2 )w^(z)^2h^ 1 [q 2 k+i(z*)c 2 k ~ q2k(z*)c 2 k+i]qi(z) 

JCr k=0 

-[q2k+i(z)c 2k - q2k(z)c 2 k+i]qi{z*) . (3.27) 
Splitting into odd and even this equation can be written in matrix form as 

A ( J = j c d 2 z{z* 2 - z 2 )w^(z) g (3.28) 

q2k+i(z*)q 2 i(z) - q 2k+ i(z)q2i(z*) q2k(z*)q 2 i(z) - q 2k (z)q2i(z*) \ / c 2 k 

q2k+i(z*)q2l+i(z) - q2k+i(z)q2l+i(z*) q2k{z*)q 2 i+i{z) - q2k(z)q2l+i(z*) J \ -c 2 k+i 

It has a solution if the following determinant of an antisymmetric 2N x 2iV matrix vanishes (after 



_ i 

normalising the q k (z) with h k 2 



= det 

l,k=0,N-l 



i 



A ^ Q fci ) - J c d 2 z{z* 2 - z>W(z) h k >h t * (3.29) 

q2k(z*)q 2 i(z) - q2k{z)q 2 i{z*) q2k+i{z*)q2l(z) - q2k+i(z)q 2 i(z*) 

Q2k{z*)q2l+i(z) - q2k(z)q2l+i{z*) q2k+i{z*)q2l+\{z) - q2k+i{z)q2l+i{z* 
N-l 

= n(*-*s fl ) • < 3 - 3 ») 

J'=0 

Here we have swapped the columns of the 2x2 matrices which is allowed under the determinant. Let 
us compare to eq. (|3.19p : using the skew orthogonality and taking out the norms which cancels the 
normalising pref actor, the matrix inside the Pfaffian there agrees with the one in eq. ([3.29)) . up to an 
overall irrelevant sign as Pf(^) = Pi{A T = -A). Because of the double de generacy of the eigenvalues 
from eq. (j3.30p we obtain as a final result for the gap probability and the generating functional 

N-l N-l 

4V) = n i 1 - a 5 4) ) . ^ (4) (^) = n - ^ 4) ) • ( 3 - 31 ) 

j=0 j=0 

We have checked that the same result can be derived independently as a quaternion determinant, using 
orthogonal quaternions instead. It would be very interesting to relate the two gap probabilities eq. 
(13.15)) for (3 = 2 and eq. (|3.31[) for f3 = 4 for a general weight function and at general non-Hermiticity. 
In the following subsection we will be able to find such a relation for maximal non-Hermiticity by 
diagonalising the determinant or Pfaffian. For the Gaussian weight we can then explicitly compute 
the A^ as a function of r. 
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3.3 Maximal non-Hermiticity (3 — 2 

In this section we derive an exact product representation for the gap probabilities for an arbitrary 
rotationally invariant weight function. The corresponding one-dimensional integrals can be solved 
explicity for our Bessel-X weight function eq. (|2.8p at fx = 1, and have been computed previously for 
the Ginibre ensemble with exponential weight eq. (|2.5f) in j5j [6] at r = 0. 

We again start with (3 = 2 and a general weight function, and consider the Gaussian chiral and 
Ginibre weights as examples at the end. For any rotationally invariant weight wo (\z\) (e.g. eq. (|2.8[) 
at n = 1) the orthogonal polynomials are given by monomial powers (we will use the chiral case with 
squared variables here, the Ginibre class trivially follows). 

This can be seen from the orthogonality relation 

/•2-ir 

/ dO z 2k z* 21 = 2irt 4k 8 kl ,k,l = 0,1,... , (3.32) 
Jo 

where z = te tS . Thus the wave functions in eq. (|3.ip only differ by their norm, defined as 

d 2 zw {2) (\z\)z 2k z* 21 = 27:5 kl / dttt ik w {2 \t) = 5 kl h k (3.33) 



with ip k (z) = h k 2 w¥\\z\)2z 2k . 

Consequently the determinant in eq. (|3.5|) becomes diagonal due to eq. (13. 32ft even though we 
only integrate over a circle of radius r. Putting the two results together we can read off the Fredholm 
eigenvalues for the gap probability 

4 2 V) = ff " *fM) . with Af (r) S • (3.34) 



j=0 



J °° dtt^W^it) ' 



After this general result let us return to our two Gaussian matrix models as examples. The weight 
function of the chiral model eq. (|2.8|) at /x = 1 contains a Bessel--fT function and all integrals can be 
done explicitly. First we compute the norms using a standard integral, see e.g. eq. (6.561.16) [20] . 



/•oo 

/ dss 2k+u+1 K u (s) = 2 2k+u {k + v)\k\ , 
Jo 



(3.35) 



after choosing the scaling variable s = N\z\ 2 . 

The second integral needed is derived in Appendix [Al 

px / 2/g I v I 1 

F v {k,x) = j o dss 2k+ ^K u {s) = 2 2k ^(k + u)\k\ (l - 22k l {k - v)m K v+1 (x) 

I [ ^\x)K u+1 {x) + I [ ^}\x)K v+2 {x))) . (3.36) 



— x 

Here we have introduced the incomplete Bessel-/ functions as 

k 



- E(T^t© 2 '"' for *- ' (3 ' 37) 

1=0 y ' 



and zero for a negative upper summation index. The normalisation eq. f|3.35|) can be written as 
F u (k,x = oo), and for the gap probability we thus have 



n 



Most explicitly we can thus write out one minus the Fredholm eigenvalues as 

1 " X k \r) = , k f v k ( l^\ K v+1 (x) + x(lf-*\x)K u+1 (x) + lf-l\x)K v+2 {x)) , (3.39) 



2 2k + v {k + v)\kV 

where x = Nr 2 . Inserted into the generating functional eq. (|3.15|) this gives all gap probabilities 
for our chiral ensemble with (3 = 2. This holds both at finite N and infinite N, where the product 
extends to infinity and the combination x = Nr 2 is kept fixed. We also note that in the large-A; limit 
the second term in eq. (|3.39p relates to the Wronsky identity for (complete) Bessel-I functions, 

= l-x(l v+1 {x)K v {x) + I v {x)K v+1 {x)\ . (3.40) 

For completeness we give the kernel at \x = 1 [18] to be inserted in eq. f|3.6[) the eigenvalues of which 
we have computed: 

Our second example is the known Ginibre ensemble at (3 = 2 [5]. Instead of eq. (|3,34p we have 

A (2) ( ) _ Jo dtt 2 ^W Gm (t) 

X > Gm[r) ~ / °° dtt 2 i +1 w Gin (t) (3 - 42) 

from the weight eq. (|2.5p at r = 0. Consequently each factor in the product for the gap probability 
can be written in terms of (upper) incomplete Gamma functions: 

1_ A ( 2 ) ( r ) - r ^ + 1 ' x ) _ e -V- (3 43) 

with x = Nr 2 . This result was first derived in [5], and the corresponding kernel in eq. (13. 6|) reads [16] 



AT Af_1 ATk 

K N (z,v*) = -e-f (N 2+ H 2 ) Y ^-z k v* k . (3.44) 

k=0 

Before turning to f3 = 4 we would like to illustrate our new result eq. (j3.38j) by plotting the individual 
eigenvalues in the chiral case, using eq. (|2.12p . For that purpose we compare the first few eigenvalues 
and their sum with the large- iV microscopic spectral density at strong non-Hermiticity. It is obtained 
from eq. fl3"ld|) [H [21] 

and has a very simple form depending only on the rescaled modulus s. For large arguments it ap- 

(2) i 

proaches the corresponding constant density of the Ginibre ensemble, PQ in {s) = -■ The corresponding 
gap probabilities are now given by an infinite product. Because of its fast convergence the individual 
eigenvalues can be computed by truncating the products at n, where we have used n = 8 in Figure Q] 
to compute the first 5 eigenvalues from eqs. (|3.16p an d (|2.1ip . Note the normalisation eq. (|2.12p of 
the Pfr(r). In Figure [T] we show tt times the density versus hpi (r). 
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Figure 1: The spectral density eq. (|3.45p of the chiral ensemble at (3 = 2 times tt, and the corresponding 
distributions of the first five eigenvalues eq. (|2.1ip . as well as their sum, for v = (left) and v = 1 
(right). 



3.4 Maximal non-Hermiticity (3 — 4 



Next we treat (3 = 4 at maximal non-Hermiticity. We choose the polynomials in eq. (|3.18|) as monomial 
powers here (these are not the skew orthogonal polynomials) as in the derivation [6] 



2 k 



d6{z* 2 -z 2 )[z 2 ^z* 2 ^ -^-i) z *2(fc-i)^ =4 vr(t 4fc - 4 5 fc _ lii -t 4fc <5 fci/ _ 1 ), (3.46) 



where z = te . Therefore the Pfafnan eq. (|3,19p can be computed as follows (see also [6] chapter 15), 
where every other term contributes, 



N-l 



3=0 



( 4 ), 



L r dt &+ 5 wi 4) (t) 



f™dt t^+ 5 wi 4) (t) 



(3.47) 



Here we have used the fact that the partition function can be computed from eq. f)3. 19|) at r = 0. 
This leads to the integrals in the denominator. Eq. (|3.47p is of the same form as eq. (|3.34p . with 
the difference being that alternate powers are skipped. This leads to the relation between the gap 
probabilities at (3 = 2 and 4 true for general rotationally invariant weight functional 



4%) 



N-l 

n 

j=0 



1 - A 



(i) 



N-l 

n 

3=0 



1 - A 



(2) 
2j+l 



(3.48) 



where only the odd Fredholm eigenvalues from (3 = 2 contribute to (3 = 4. The same statement can be 
made relating the generating functionals, as well as for the Ginibre ensembles. This new relation can 
be compared to the relation known for Hermitian Gaussian RMT (r = 1) in the large N limit where 



one has E^\r) = ±( TlgoU " + II£o(l " A 



(2), 



,(2) ' 
2i+l. 



given as the sum of the product of even 



and odd eigenvalues (see e.g. [6] eq. (11.7.5)). It would be very interesting to generalise our relation 
eq. (|3.48p to the general case where one could take limits /i -> or [i — > 1. We can only conjecture 
that ^(r)-dependent prefactors of a general linear combination of the product of even and odd terms 



We have to use the same weight for both [3 here, that is keeping [3v fixed in eq. (|2.8 
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could provide a form valid in the limiting cases. For completeness we also give the expression for the 
Fredholm eigenvalues in the Ginibre ensembles, 



A 



(4) 

j Gin 



(r) 



ft; dttv+ 3 w Gin (t) 

/ °° dtt 4 i+ 3 w G in(t) ■ 



(3.49) 



Switching back to our explicit example of the chiral Gaussian model with a Bessel-.fr weight we 
do not need to do a new computation and can simply use the integrals provided by eq. (I3.39P 



N-l 



n i 



fc=0 



F 2u (2k + l,Nr 2 
F 2v (2k + l,oo) 



(3.50) 



In comparison to (3 = 2 we have to shift the index v — ► 2v here, due to the explicit /3-dependence of 
the weight eq. (|2.8p . The corresponding pre-kernel to be inserted into eq. (|3.23j) is given by [8] 



K N (Z,V*) 



N 2u+2 N-l k]T ^ + y + l)N 2k+2j+l ( z Ak+2 v *ij _ z 4^*4fc+2) 



o p^ Q T(2k + 2u + 2) (2k + 1)! 2^jT(j + v + 1) 



(3.51) 



We also give the corresponding expression for the Ginibre ensemble at f3 = 4 which was already 
derived in [6], 



1 " >fUr) 



1 - A 



(2) 

2j+l Gin 



(r) 



T(2j + 2,x) 
T(2j + 2) 



E 

k=0 



X 



with x = Nr 2 . The pre-kernel is here given by [T7] 



AT— N-l k 

N 2 



2vr ^ ^ (2fc + l)!!(2Z)!! 

fc=0 z=o v 7 v ; 



2/c+l„,*2Z „2Z„,*2fc+l 



V 



Z V 



(3.52) 



(3.53) 



S? P f\se*)de 




Figure 2: The integrated spectral density eq. (|3.54p of the chiral ensemble at (3 = 4, and distributions 
of the first five eigenvalues eq. (|2.1ip . as well as their sum, for v = (left) and v = 1 (right). 

As before we illustrate the new individual eigenvalue distributions we found for (3 = 4 in the chiral 
case resulting from eqs. (I3.16|) and (I2.11|) . The microscopic spectral density in the large- N limit at 
maximal non-Hermiticity is given by [8] (note the radius of support (3/2 and height depend on (3) 



^)H 2 K 2v (M 2 ) J 1 -jB= I 2! ^M 2 )sinh (hjT^&tf - rf 2 )) , (3.54) 
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where 77 = \^Nre ld . In contrast to /3 = 2 the density is no longer rotationally invariant and explicitly 
depends on the angle 9. In particular it vanishes along the real and imaginary axis, and for a more 
detailed discussion of its symmetries we refer to [8]. In order to be able to compare with the individual 
complex eigenvalues we have to integrate over the angle 9 as indicated in eq. (|2.12p . The result for 
the first few eigenvalues is shown in Figure [2 Here we have again truncated the infinite product at 
n = 8 to display the first 5 eigenvalues and their sum. 

4 Asymptotic expansion for the gap probability for Ginibre (3 = 2 

(2) 

We now turn to consider the large-r asymptotics for the gap probability Eq (r) in the large N limit 
keeping x = Nr 2 fixed. Our method has the advantage that it can be easily extended to the chiral 
ensembles in the next section, although for this part it provides a weaker result than given in |13j . 

Here we use eq. (j3.34j) , but with the Ginibre Fredholm eigenvalues from eq. (|3.43p . Starting at unity, 

(2) 

as r increases E (r) tends towards zero. Therefore, we define P(x) > by 

Ef\r) = exp [-P(x = Nr 2 )} , (4.1) 

and will analyse the large- x behaviour of P(x), which is given by 

00 00 

P(x) = -loglJ^O*) = -Y, l °S9n(x) (4.2) 

n=0 n=0 

where 

n y, 

9n(x) = ^2t k (x) and t k (x) = e~ x jj- (4-3) 

k=0 

The large N limit of the corresponding Fredholm equation f|3.6|) reads 

\^(rj) = r dtt d9 e-^ t2+ ^ 2)+r)te ~ ie ^{te ie ) . (4.4) 
Jo Jo 

4.1 Statement of large x asymptotic results 

We show that using our method 

N X 2 xlogX /fog27T \ , — 

P(x) = — + ^^ + l^--l)x + 0(Vx) (4.5) 

although the coefficient of x is partly conjectured. It was derived rigorously in [13] including the co- 
efficient of the next order term. In more detail, we proceed by splitting the sum over n in eq. (|4.2j) 
into two parts, namely S\(M,x) (in which n runs from 1 to J\f) and £2 (AT, x) (with n running from 
M + 1 to 00), where J\f is some arbitrary (at this stage) parameter (which in general can depend on x) 
0. We then further split S\{N, x) into two parts, denoted Su(J\f,x) and S±2 (A/", x), so that we have 

P(x) = S 11 (M,x) + S 12 (M,x) + S 2 (N,x) 

MM 00 
= - ^ log t n (x) - ^2 lQ g r n (g) ~ ^2 l °S9n{x) (4.6) 

n=0 n=0 n=M+l 

4 7V here is unrelated to the matrix size JV which has been taken to infinity. 
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with r n (x) defined as 



. . Qrt(x) n n(n — 1) n\ 

r n {x) = = 1 + - + 1 „ ' +- + — 



t n {x) 



X 



(4.7) 



and t n (x) and g n (x) denned above in eq. (14, 3ft . Since Si2(M,x) < 0, we will always work with 
\Si2(M,x)\ in the following. 

With the choice that M = x (which we justify in Subsection 14. 5h . we show that the following three 
properties i) - in) hold: 



i) Sn(x,x) = — + 



X 2 X log X log 2"7T 



+ x + 0{\ogx) 



(4.8) 



ii) 5^(x,x) < |5i 2 (x,x)| < Stf(x,x) 



UB, 



where 



and 



SY 2 B (x,x) = x + 0(\ogx) 



$12 ( X 7 x ) 



l-lo g 2 + ^ + iLi 2 (-e- 2 ) 



with Li 2 (x) being the dilogarithm 



OO I. 



Li 2 (x) = ^^2 



fc=l 



(4.9) 



(4.10) 



z + 0(l) « 0.653a; + 0(1) (4.11) 



(4.12) 



) S 2 (x,x) = O(y^) • 



(4.13) 



More precisely, 
where 

and M is a constant given by 



< S 2 (x,x) < S% B (x) 
S% B {x) = My/x + 0(l) 



M = - I log(*(m))dm w 0.478 
'o 



and $(m) is the cumulative normal function 



1 f r 



(4.14) 
(4.15) 

(4.16) 
(4.17) 
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Furthermore, based on some preliminary numerical analysis, we conjecture the following somewhat 
stronger result for S12 (we comment on the bound for the sub-leading term later) 





ii)' -Su{x,x) = X + 0(y/x) 


(4.18) 


and for S2 we conjecture 




in)' S 2 (x,x) = My/x + C + o(l) 


(4.19) 



where M is defined as before (eq. (|4.16p ). and C is another constant, given by 

1 r°° (m 2 - l)e~ m2 / 2 



C = 



-dm — log 2 



GV^tt Jo $M 
4.2 Calculation of Sn 

We begin by computing Su containing the first two leading orders, 



-0.716 . 



(4.20) 



Sn(M,x) 



n=0 



g e 



71! 



^^(x — nlogx + logn!) 



n=0 



AT 



... Af (Af + 1) . v^, 
(A/ + l)x ^— log x + log n! 



(4.21) 



n=l 



(since logO! = log 1 = 0). We evaluate the sum of factorials as follows: 

J\f J\f n J\f M 

^logn! = E^logfc = ^(AA - k + 1) log k = (M + 1) £ log n - £ n log 



AT 



n 



n=l 



n=l k=l 



k=l 



n=l 



,i=l 



(M + 1) log ATI - n log n. 



n=l 



For the first term, we use Stirling's formula 



and for the second term, we apply the Euler-MacLaurin summation formula (as hinted in [5]) 

N 

£/(») = 



n=l 



(4.22) 



(4.23) 



(4.24) 



which gives 
y^ n log n 



n=l 



tlogtdt + " 8 + O(logAA) 



1 2 
AA 2 logAA AA 2 AfloeM 



r t 2 logt t 2 
2 7 



A/' 



+ ^— + 0(logA/) 



4 + 2 



+ 0(logAA). 
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(4.25) 



Inserting eqs. (14.23|) and (|4.25j) into eq. (14.22|) gives us 

1 

3M 2 ... /log2vr 



S u (AF,x) = (Af+l)x -logxH — 



+ AHogA/ r + —2- 1 AT+0(logAT) (4.26) 



4 

and on setting M = x, we arrive at 

r, / \ xlogx log27r . . . 

5ii(s,z) = — + ^^ + ^— x + 0(logx). (4.27) 

4.3 Calculation of 5i2 

We turn now to the large x behaviour of 

M 

\S 12 (M,x)\ = Y,^gr n (x), (4.28) 

71=1 

with r n (x) defined in eq. (|4.7p . We can write the sum starting at n = 1, because logro(:c) = 0. 

Unfortunately, we have found no simple method for calculating S±2, and so we determine some 
upper and lower bounds, and then conjecture the behaviour of S12 based on some numerical investi- 
gations. Note that, in our notation, (^19^) is <^ upper (lower) bound for 15*121 (rather than for 
S12 itself), since S12 < 0. 

4.3.1 Upper bound for | *S'i2 1 

Let us first obtain an upper bound for \Si2(M,x)\, assuming N < x. We have 

, , n n(n — 1) n\ 

r n (x) = 1 + 2 +■■• + — 

n rn\ 2 /n\ n 1 

< 1 + ^ + y + +•■■ = —- ( 4 - 29) 

X 

for n < x. Therefore 

logr„(x) < -log(l-^) (4.30) 
from which it follows that (for M < x) 



\S 12 (M,x)\ < S$(M,x) = - J^log (l - ^) . (4.31) 

n=l 



We approximate the sum using the Euler-MacLaurin summation formula. Since our analysis is valid 
only for J\f < x, we choose M = x — 1, and this gives 

Sff(x-l,x) = x-^logx + 0(l). (4.32) 

We consider the term n = N{= x) separately; each term in the definition of r n {x) is individually < 1, 
and there are n + 1 such terms, so 

< logr n (x) < log(n + 1) = O(logx). (4.33) 

Hence we have 

S${x,x) = x + 0(\ogx). (4.34) 
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4.3.2 Lower bound for |5i 2 | 



We have not been able to find a particularly tight lower bound for |<Si 2 | using elementary methods. 
Our best attempt involves the use of the following lower bound for r n {x) (denoted r^{x)) which is 
strong for low values of n (indeed, we have equality for n = 1 and 2), but rather weaker when n is 
closer to (but still less than) x: 



, . n n(n — 1) n\ 

r n (x) = 1 + „ ' +...+ 



it 



x 

„LB 



rf{x). 



Therefore 



logr-^(x) = log I 1+11 + 



2\ n 



log 2 = log 1 + 



+ log 1 + 1 + 



n 



>o g (^VE H ' +1 



k=l 



Hence 



k 



2 

1 + - 

x 



-kn 



log 2. 



(4.35) 

log 2 
(4.36) 



Si 2 (A/» > Sg(AT,x) = ^logr^(x) 



n=l 



M(M + 1) 



M oo 



(-) 



fc+1 



k 



1 + 



»=1 k=l 

It is convenient at this stage to set M equal to x. Using the result that 



-kn 



AHog2. (4.37) 



1 

1 + - 

a 



e + 0[ - 

a 



we can easily evaluate the first term of eq. (|4.37p 



x(x + 1) 



log 1 + 



(x + 1) log 1 + 



x/2 



(x + 1)(1 + 0(x~ 1 )) 



= x + 0(l). 

The second term, which can be evaluated using the result (which we do not prove here) 

— kn 

N J 



N oo 

EE 

71=1 fc = l 



( _ )fc+ i f 2_y 



k 



V -Li 2 (-e 2 ) 

24 2 1 ' 



is also linear, as indeed is the third term (which is simply — x log 2) 
Combining these results then gives 



"^12 ( x i x ) 



l-log2 + ^ + iLi 2 (-e- 2 ) 



x + 0(l) 



(4.38) 



(4.39) 



M + 0{1) » 0.346AA + O(l), (4.40) 



(4.41) 



where Li is the dilogarithm defined in eq. (|4.12p . The term in square brackets is approximately equal 
to 0.653. 
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4.3.3 Conjecture for |5i2| based on numerical analysis 

In the absence of a concrete analytical proof for the asymptotic limit of S12 itself, we undertook some 
elementary numerical analysis. In particular, we considered 

x 

for various (necessarily finite) values of x. The function A\2{%) begins at ~ 0.69, and is increasing. 
It flattens very quickly, and appears to be converging monotonically to 1. For example, we find that 
yli 2 (200 000) m 0.996. 

Assuming that the coefficient of x is indeed unity, we then numerically investigated 

s \S 12 (x,x)\-x^ 



for increasing x. It was not completely apparent that this converges (we again went as far as x = 
200 000), but it seems possible that it does, implying that the sub-leading term is 0(y/x). 
Putting these together, we anticipate that 

\S 12 (x,x)\ = x + 0{y/x), (4.44) 

(i.e. that the true asymptote actually equals our earlier upper bound, at least to leading order). 
Considering the sub- leading term, this might (at first sight) appear inconsistent with eq. (|4.34[) . since 
we have replaced O(logx) with 0(y/x), but there is no contradiction. The more complete statement 
of our conjecture is, in fact, that 

\Su(x,x)\ = X — C\[x + o(y/x) (4.45) 

where c is some positive constant. 

4.4 Calculation of S 2 using an upper bound 

We need only to provide an upper bound for £2(2;), since our aim is merely to show that S^sc) is 
smaller than linear in x, and we know that S 2 (x) > 0. 
From eq. (|4.6|) we have 

00 

S 2 {M,x) = - ^gg n (x). (4.46) 

n=Af+l 

For n > x, we have (see Appendix [B] for details) a lower bound for g n (x): 



g n (x) > g ^( m ,x) = ^m)-J2^W--^e-^-^^ (4.47) 

. -, X ' V X 

1=1 

where we have introduced a scaled variable 



n — x , . 

m = — t=- (4.48) 
\Jx 

and <I>(m) is the cumulative normal function. The Ci(m) are m-dependent numbers, and a, (3, 7 and 
S are constants. Therefore we have (on setting N = x) 

S 2 (x,x) < S^ B (x) 

00 />oo 

log(5 LB (m,x)) dn = -y/x / log(g w (m,x)) dm (4.49) 
x Jo 
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and on factorising out the $(m) and then expanding the logarithm, we arrive at 

/•oo 

S% B (x) = -y/x log($(m)) dm + 0(1) = My/x + 0(l) (4.50) 

J 



where 

M 



POO 

- / log($(m)) dm ps 0.478 , (4.51) 
Jo 

see also eq. (27) in [13] . For our purposes, it is sufficient that we have shown S2(x,x) = 0(y/x). 
However, some more detailed analysis (supported by numerical analysis) leads us to conjecture the 
following, much stronger, result: 

S 2 {x,x) = MV5 + C + o(l) (4.52) 

where 

1 f°° (m 2 - De'™ 2 / 2 

C = —= / ^ — ^ dm -\og2 » -0.716. (4.53) 

6^ Jo $(m) 

4.5 Choice of the 'split point' J\f 

Our final result must, of course, be independent of the choice of J\f. However, if we had chosen N <C x, 
then we would not have been able to bound S2 in the way that we did, and it would not have been 
of smaller order than Si. Conversely, if we had chosen M 3> x, then our argument for S12 would not 
have been applicable, and this term would then have been of higher order. So, setting N = x is a 
pragmatic choice which simplifies our analysis, giving Sn S> IS12I 

We can understand this in a little more detail if we analyse what happens to each of Sn, S12 
and S 2 as we increase the value of N by 1. For Su, we can estimate the effect of this by partially 
differentiating eq. (I4.26|) with respect to M: 

aa 1 

—± = N (logAT - logx) - - x) + logAT - -logx + 0(1). (4.54) 

We can effectively 'minimise' this by setting M = x, which has the effect of killing everything down to 
O(logx). So, by setting M = x + 1 instead of M = x, we will only see an increase in £11 of O(logx), 
which we have already decided to call 'small'. 

Of course, S12 must decrease by O(logx) to compensate, and the argument preceeding eq. f|4. 33|) 
shows that this is precisely the case. (S2 will change only by 0(1), and so can be ignored.) 

So, with this choice of M, we can be sure that all of the 'large' contributions are already included 
in Su; consequently, it is the optimal choice. 

5 Asymptotic expansion of Eq(t) for the other ensembles 

We now extend the problem to the j3 = 4 case, and to the corresponding chiral ensembles. For 
completeness, we will also include here the (3 = 2 Ginibre case which was discussed in fuller detail in 
the previous section. 

The aim here, as earlier, is to determine P(x) where E${r) = ex.p(—P(Nr 2 )). We established 
before that 

n 
fc=0 



X 



2n+v+l Ku+l{x) 



< n{x) = ^nl(nl7)\ +'[Ku + i(x)I^ 2 %) + K v+2 {x)I [ ?- 1 1 \x)\ (5.2) 
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Ensemble 


(3 = 2 


(3 = A 


Ginibre 


oo 

P{x) = -Y J ^ in {x) 

n=0 


oo 

P(x) = - Yl ^£ in {x) 

n=l, n odd 


Chiral 


oo 

P(x) = -J>g 

n=0 


oo 

P(X) = - Y lQ g<n(^) 
n=l, n odd 



Table 1: Definition of the exponent P(x) for all four ensembles. 

where -Zj^(x) is the incomplete /-Bessel function defined in eq. (I3.37|) . 

Our task is to find the asymptotic behaviour of each of the P(x) which have the form given in 
Table [H For (3 = 2 chiral the large N limit of the Fredholm equation (13. 6p takes the compact form 



2k 



M{V)= dtt dd t\ V \K u {t z )* K„(\ V \ Z )2 I u { V te- l °)^(te w ) , 
Jo Jo 



(5.3) 



after inserting the radial Bessel kernel. For (3 = 4 a similar equation can be written, but the corre- 
sponding kernels from eqs. (|3.5ip and (|3.53j) do not simplify. 



5.1 Statement of results 

To proceed, we first identify the leading term (highest power of x) of each g n {x) 

tf n {x) 



111 



{ x\2n+v+^ 
-x\2> 2 



n\(n + u)\ 



and denote the 'remainder' as 

„Gini 



r u,ni x ) 



9u h n( x ) 



n n(n — 1) n\ 



We have used the asymptotic form for the iT-Bessel function 



K u (x) - J— e 



(5.4) 
(5.5) 

(5.6) 
(5.7) 

(5.8) 



as x — ► oo when determining tf^ n (x) in eq. (|5.5|) ; however, the terms that we dropped in doing this are 
incorporated into r„ n (x), so no approximation has been made. For notational convenience, we drop 
the l Gin/ch' 1 labels in what follows, and write everything generically. 



22 



Ensemble 


Choice of 


Sn(M,x) 




Coeff. of x in P(x) 


Ginibre (/3 = 2) 


X 




A(x) 


— X 


log 2vr 
2 


Ginibre (/3 = 4) 


X 




A(x) x 
2 4 


X 

~ 2 


log 2vr 3 
4 4 


Chiral (/3 = 2) 


X 

2 




A(x) v 

X 

2 2 


X 

~ 2 


log 2tt v + 1 
4 2 


Chiral (/3 = 4) 


X 

2 


A(x 
4 


) (2u 1\ 


X 

~ 4 


log 2vr 2^ + 2 
8 4 



Table 2: Explicit results for the exponent P{x) for all four ensembles. 



We can write P{x) as the sum of three sums: 

Sii(A/» = - J] lo §^( x )> ( 5 - 9 ) 

0<n<A/" 

Sn(N,x) = - l °Sr n (x), (5.10) 

0<n<A/" 

S 2 (Af,x) = -J2 l °S9n(x) (5.11) 

where N is a suitably chosen 'split point', and n is restricted to the odd integers for the (3 = 4 
ensembles. To simplify the appearance of our results, we introduce A(x), defined as 

x 2 xlogx log2vr 

Our results are then summarised in Table [2] where we omit terms of smaller order than linear in x. 
For Sii these terms are O(logx), and for S12 they are 0(^/x). In all cases, S2(x) = 0{yfx). 

5.2 Derivation by relating the ensembles 

The method for each of the 3 other ensembles is broadly similar to that for the original (Ginibre (3 = 2) 
case, and so we do not give details here. Essentially, the same techniques can be used to determine 
the relevant bounds, and we have undertaken numerical analysis to support the results that we claim 
above. 

However, it is quite instructive to understand the relationships between the results for the different 
symmetry classes, so we will provide here some alternative (and quite concise) proofs of the S11 's 
which highlight these relationships. 
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5.2.1 Prom Ginibre to chiral 

Assuming the Ginibre results, we will show how to derive the corresponding chiral results (for v = 
only), starting with the (3 = 2 case. We begin by stating the following exact relationship which follows 
immediately from the definitions: 

1 , x 1 



Therefore 



log t* (x) = 2 log tf n (|) + - log | + - log tt. (5.13) 



S^=\ X ) EE - kg *?»(*) 

0<n<rr/2 

= "2 E log^ff)-^ £ (kg f + log- 



25: 



2/ 2 ^— ' V 2 

0<n<x/2 0<n<x/2 
Gin,/3=2 ( x 



2 £ ( log 1 + log 7r 



11 V27 2 

0<n<x/2 

= 2A(|)-|(log|+log7r)+0(logx) (5.14) 

where A(x) was defined in eq. (I5.12|) . But it easily follows from the definition of A(x) that 

2A fx\ = M^ + x / x + log? x + 0(lQgx) (5 _ 15) 



2 4 V 2 

where the scaling of the leading term is explained by the fact that A{x) is essentially quadratic in x. 
Hence we have 

5,r= 2 (x) = ^ + 0(logx) 

qGin,f3=2 / \ 

= 11 2 {X) +0(logx), (5.16) 

i.e. the chiral sum is half the Ginibre sum. This is the result given in Tabled (when v = 0). 
A similar argument gives the corresponding result for (3 = 4 (at v = 0): 

qGin,f3=4:/ \ 

S^=\x) = 11 2 (X) -|+0(logx). (5.17) 

Note that here we pick up an extra term which is linear in x, so it is not true in general that the chiral 
case is always half the Ginibre case (at least, not when we consider the terms linear in x). 



5.2.2 From (3 = 2 to p = 4 

The (3 = 4 case involves summing alternate (odd n) terms from the same sequence used for (3 = 2. We 
would therefore expect the total for (3 = 4 to be approximately half that for (3 = 2, since — log g n (x) is 
a 'smooth' function of n (for fixed x). However, since — log <? n (x) is monotonic decreasing as a function 
of n, we would expect the sum of the odd terms (starting at 1) to be slightly less than the sum of the 
even terms (which start at zero), and indeed we do see this bias. 
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Let us now quantify this; we will do this first (and in detail) for the Ginibre ensembles. We have 

0<n<x ^ ' ' 

= - Yl {(y loge ~ X J^T[fi) + log ( e ^)} ( + if x even) 

n odd 

S { 2 ( l0ge ^^) " l0g (3} < + ditt °) ( 5 - 1§ ) 



l<n<x 
n odd 



and thus 



S f^= 2 (x) = 2S^=\x) + ]T log- (+ ditto) 



l<n<x 
n odd 



2sfi n >P- 4 ( x ) + Y ^gx- Y log(2n) + 0(logx) 

0<n<x/2 0<n<x/2 

2S^ =A {x) + ~ +0(logx) (5.19) 



where the last summation has been approximated using the Euler-MacLaurin summation formula. A 
trivial rearrangement then gives^l 

qGin,f3=2 1 \ 

Sfr" = %) = 11 2 ( " j - \ + OQogx) (5.20) 
which is the result given in Tabled 



An almost identical argument gives a similar result for the chiral ensembles (for v = 0) 



5. 



5^ = 4 (x) = ^ - (X) - | + 0(log x) . (5.21) 
5.3 Illustration of our results 

Writing eq. (15. ip in terms of gamma functions (see eq. ()3.43p ) allows us to generalise g^ m (x) to 
non-integer n. For the case when x is large, we can then easily show that 

fyx) ^ g^ +v+k {x) . (5.22) 

The P{x) for all of the four ensembles can therefore be written as products of the logarithms of certain 
g^ m (x) as given in Tabled 

We can easily see from this why the Ginibre (3 = 4 case is slightly less than half of the Ginibre 
(3 = 2 case, for example. There are half as many dots in the second row as in the first, but they are 
systematically shifted to the right (this is the same argument we presented above in Section [5.2.2p . 

More interesting is the comparison of the chiral and Ginibre cases for fixed (3. For the chiral [3 = 2 
case (for v = 0, i.e. no exact zero eigenvalues) there are also only half as many dots, but there is no 
systematic shift compared with the Ginibre (3 = 2 case as the chiral dots take the average value, e.g. 
2 = (0 + 1/2) etc, and so the chiral S\{x) is exactly half that of the Ginibre case. 



5 We could have changed the proportionality factor from | — > 1 by modifying the weight to exp[— ^iVTr . 
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Ensemble 




4 - 5 — 

* 2 ° 2 


1 Q 

6^7 etc. 


Ginibre (3 = 2 


• • • • 


• • 


• • 


Ginibre (3 = 4 


• • 


• 


• 


Chiral (3 = 2 (v = 0) 


• • 


• 


• 


(u = l) 


• • 


• 




(u = 2) 


• 


• 


• 


Chiral (3 = 4 {v = 0) 


• 




• 



Table 3: Pictorial illustration of the Fredholm eigenvalues taken by the four ensembles. 

What about the case where we know there is exactly one eigenvalue at the origin? For the Ginibre 
(3 = 2 case, we must then consider the probability that, given there is one eigenvalue at the origin, 
there are no other eigenvalues within a distance r of the origin. It turns out that this conditional 
probability is calculated simply by removing the first (n = 0) term from the sum that defines Si (see 
also [5] ) . Compare this modification of the top row of Table [3] with the chiral (3 = 2 case for v = 1 
(i.e. when there is again precisely one exactly- zero eigenvalue). It will be seen that there are half as 
many dots in the latter case as in the former, and there is again no systematic shift. So, for precisely 
one exact zero eigenvalue, it is also the case that the chiral case is exactly half the Ginibre case. 

It is not possible to do a similar comparison for (3 = 4 directly, since the Ginibre ensemble has zero 
probability to find an eigenvalue at the origin. 

5.4 Low x asymptotics for the four ensembles 

The probability that one eigenvalue is at the origin and a second one at radial distance r can be 
used to compute the spacing distribution p(s) in the complex plane [5]. Following their argument the 
mean of the macroscopic large- N density is constant on a disc for the Ginibre ensemble (3 = 2 and 
one can assume that the probability calculated in this way is translation invariant, giving the spacing 
everywhere in the bulk. 

However, for all the other ensembles this is not true and the origin is special. What can be 
compared is the strength of repulsion between two complex eigenvalues to leading order: for the 
Ginibre ensemble (3 = 2 everywhere in the bulk following [5] , and for the chiral ensembles at the origin 
by placing an exact zero eigenvalue at zero. 

Therefore we also give the expansion of the gap probabilities close to the origin. Note that the 



Ensemble 


(3 = 2 


(3 = 4 


Ginibre 


x 3 5x 4 7x 5 

1 - x H 1 h ... 

2 12 24 


2 S 4 ^ 
kJC *Xj tJC CjC> 

1 1 1 h ... 

2 3 6 15 


Chiral [y = 0) 


1 + — log x + . . . 


1 + T77 log X -\- . . . 

lb 


(u = 2/(3) 




4 

X 

1 "64 + - 



Table 4: Expansion of the gap probabilities for small radii x = Nr 2 <C 1. 
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first product begins at zero, unlike in [5], as we give the gap probability here in Table |4j compared 
to the conditional probability to have one eigenvalue at zero and one at radius r there. The first two 
equations agree with [6] eqs. (15.1.19) and (15.2.18). 

The leading power of the spacing distribution p(s) can be reproduced for Ginibre = 2 as follows. 
Multiplying the first line in Table U] with e x to remove the n = contribution also removes the linear 
term — x, and we obtain 1 — ^x 2 + . . .. The spacing is then obtained by first reinserting x = Nr 2 and 
then differentiating with respect to r: p(s) ~ s 3 . The same cubic repulsion was found in [22] for other 
examples of rotationally invariant weights, as well as in [16] for strong non-Hermiticity by expanding 
the gap to second order as in eq. (|2.15p . and it is thus considered to be universal. Comparing this to 
the chiral ensemble with one exact zero eigenvalue v = 1, we find once more a cubic repulsion, adding 
a further ensemble to this universality class. 

While for (3 = 4 a comparison to Ginibre is not possible we can at least compare the two chiral 
ensembles, finding that the repulsion is much stronger for (3 = 4 than for (3 = 2: here placing one 
eigenvalue at the origin corresponds to v = h. This indicates that there exists a different universality 
class p(s) ~ s 7 of spacing distributions for (3 = 4 at the origin. 

As a final observation we note that when comparing ensembles in Table H] the first power in x 
doubles when going from Ginibre to chiral, or from /3 = 2 to (3 = 4. This is in contrast to the coefficient 
in the exponent going down by ^ for these comparisons of ensembles in the large- x asymptotics in 
Table [2 

6 Conclusions 

We have investigated the gap probabilities and distributions of individual eigenvalues with respect to 
radial ordering in the complex plane in non-Hermitian Random Matrix Theory (RMT). 

After setting up a general framework in terms of Fredholm determinants and Pfaffians for general 
non-Gaussian RMT with unitary {(3 = 2) and symplectic {(3 = 4) invariance we turned to maximal non- 
Hermiticity. For general weights with rotational invariance we found that the product representation 
of the gap probabilities in terms of Fredholm eigenvalues are related for (3 = 2 and 4, and we gave 
explicit new expressions for two ensembles of Gaussian chiral RMT. This relation between the gap 
probabilities is different from the one for Hermitian RMT. It would be very interesting to extend 
our relation to intermediate Hermiticity, interpolating between the two limiting cases. This may be 
compared to the Hermitian limit of the spectral correlators. For (3 = 4 these are given in terms of 
a 2 x 2 matrix kernel containing a single pre-kernel of skew orthogonal polynomials in the complex 
plane which is much simpler than the three pre-kernels in the Hermitian limit following from a Taylor 
expansion of the complex one. 

We then derived an asymptotic expansion for the gap probability at large radii for the Gaussian 
Ginibre and chiral complex ensembles, both at (3 = 2 and 4. In particular this included a detailed 
discussion how to get from Ginibre to chiral and from (3 = 2 to 4 in these ensembles. Our results 
are consistent with the known results for Ginibre at (3 = 2 of Forrester, for the other ensembles they 
were new. It would be very interesting to sharpen our strict upper and lower bounds for the linear 
coefficients, and we have given numerical evidence for their conjectured values. Expanding for small 
radii we found that the chiral complex (3 = 2 ensemble also displays a cubic level repulsion at the 
origin, in contrast to its (3 = 4 counterpart. 

While in this paper we have focused on the (3 = 2 and 4 ensembles one could try to generalise 
our results to the recently solved (3 = 1 Ginibre and chiral ensemble with real asymmetric matrix 
entries. Because these two ensembles have both real and complex eigenvalues, several different gap 
probabilities can be defined and computed, and we leave these open questions for future work. 
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A Integrals over Bessel functions 

In this appendix we derive several integrals over Bessel functions, including the one given in eq. (|3.36p . 



A.l General case 



(21 

First we compute the matrix elements of a determinant yielding the gap probability E (r) in the 
general elliptic case of the chiral ensembles. For that purpose a slightly different representation than 
eq. (|3.ip is more convenient. Instead of choosing orthonormal wave functions in eq. (|3.ip one can 
simply keep the original monomials in the Vandermonde determinant, leading to 



m 



det 



cr( 2 ) l<k,j<N 
ch 



C\Cr 



d 2 zw^(z)z 2 ^ k - 1 h* 2 ^ 



(A.l) 



instead of eq. f)3.5j) (the same identity was shown in a different way in [6] eq. (15.1.13), without squared 
arguments). While this seems tailored to the rotationally invariant case at maximal non-Hermiticity 
[U] we can also use it here in the general case. We have to compute the following integrals: 



2 77 



dss I d6s 2u+2 K u (as 2 )exp bs 2 cos(26) 



poo 

7T / dU fe+ ^- 1 ^(ai)/ fc+i _ 2 (6i) 

Jr 2 

N(l + n 2 ) , N(l-n 2 ) 
with a = \ S 1 , b = 



a 2(fe+j-2) 



exp[2i9(k + j - 2)} 



2^ 



2^ 



(A.2) 
(A.3) 



and we have a > b > ensuring convergence. While we were unable to give a general result we can 
provide a recursive prescription, starting with the known integral 



dttK (at)I (bt) 



1 



i 2 -& 2 



fer 2 Ii(6r 2 ) J PC (ar 2 ) + ar 2 /o(6r 2 )Ki(ar 2 ) = f(a,b;r 2 ) , (A.4) 



for v = and m = k+j — 2 = 0. Any value of m > and v > can be obtained from this by applying 
the following Bessel identities for a, b ^ 



m 



d b I m (bt) - —I m (bt) 



-d a K v {bt) + -K v {at) 
a 

We thus obtain as a final result for m = k + j — 2 

/•oo 

/ dtt m+u+1 K u (at)I m (bt) = 

Jr 2 



tl m+1 (bt) 

tK u+ i{at) . 



(A.5) 



On 



m 



m — 1 



d b -d a + 



d a + 



v - 1 



..(-d a )f(a,b;r 2 ) . (A.6) 
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All differentiations can be carried out recursively and are algebraic, acting on the expression in the 
middle of eq. (|A.4p . For b = or \i = 1 we have rotational invariance and all integrals can be 
computed in a closed form, see Subsection 13.31 and this appendix below. 

We end this part on the elliptic case by sketching how to proceed for (3 = 4. Keeping monic powers 
instead of skew orthogonal polynomials in eq. (|3.19p we arrive at 



^0 \ r ) - jr( 4 ) ^h<k,l<2N 



j2 I 2 *2\ (41/ \( 2k-2 *2l-2 *2k-2 2l-2\ 

a z[z — z jwl'yzjiz z —z z J 

C\Cr V J 



(A.7) 



After multiplying out we have to compute the same types of integrals as already done for = 2, with 
the index 2v of the Bessel-.fr function and weight now increasing in steps of two instead. 

A. 2 Maximal non-Hermiticity 

In the second part of this appendix we compute the integral needed in the rotationally invariant case 
fi = 1 (b = 0) given in eq. (|3.36[) . In this case the angular integration diagonalises the matrix, see 
(jA.ip . and the Bessel-I functions from above become simple powers. We have to show that: 

F v (k,x) = f dss 2k+y+l K v (s) (A.8) 

„2k+v+l 







2 2k ^{k + u)\k\ ( 1 - 2 2kL k + „ VM K » +l{x) (A - 9) 



k-2 „ „, , , „ k-1 



where the sums ^i=qi Tld=o are se ^ to zero - For k = this integral is standard, see e.g. eq. (6.561. 
in \M 



[ X dss u+1 K u (s) = 2 y v\-x u+1 K u+l {x) , (A.10) 
Jo 

and in this case the sums giving the incomplete Bessel-I functions in the second line are absent. Using 
the identities 

f> +1 /^ +1 ( S ))' = -s u+l K u (s) and (s^K^s))' = - S - V K V+X {s) , (A.ll) 
we can show that the following recursion holds: 

F„(k + l,x) = - f dss 2k+2 (s u+l K u+1 (s))' (A.12) 



Jo ' 
= -x 2k+2+u+1 K u+1 (x)-2(k + l) J X dss 2k+2+2u (s- v K v {s))' 

= -x 2k+3+ »K u+l (x) - 2{k + l)x 2k+2+ »K v {x) + A{k + l)(k + 1 + v)F v (k,x) . 
It is easy to verify that the explicit expression eq. (|A.9P satisfies this relation 



2 2k ^(k + 1 + „)!(* + 1)1 ( 1 - 2 n +2+u * k + 1 + v)Kk + 1 y K v+1 (x) 



2k+2+u+l 



1=0 
k-1 



1 /gx «+«/+! j fx ^2k+u+l 

- X 2^ {l + v + 1)m UJ K » +2{x) - X (k + v + l)!k! uj K ^ 
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In a last step we have to apply the Bessel-if identity 

2(i/ + 1)K V+1 {x)-xK v+2 {x) = -xK„(x) , (A.13) 

to arrive at 

F u (k + l,x) = A{k + l + v){k + l)(F u {k 1 x)+x 2k+v+1 K u+l {x)) -x 2k+3+u K u+1 (x) 

-A(k + l)kx 2k+u+1 K u+1 (x) - 2(k + l)x 2k+u+1 K v+2 (x) (A.14) 
which finishes our proof by induction. As before the same integrals computed here apply to f3 = 4. 

B A lower bound for g^ m (x) when n > x 

We derive a lower bound for g^ m (x), valid when n > x. In order to do this, we first need to determine 
an upper bound for the individual terms in g^ m (x), denoted tfc(x), for k > n > x. We will actually 
do this for k > x. 

B.l An upper bound for the reciprocal of a factorial 

We have that [23] 



n\>V2Kn^-J (B.l) 

1 



and thus 

B.2 An upper bound for tk(x) 
We define 

tk(x) = (B.3) 

and wish to determine an upper bound for this when k > x. It proves useful in what follows to 
introduce a scaled variable 

k oc , . 

P = —?=- > 0. (B.4) 



'x 

We use eq. (1B.2|) . and write k in terms of p as follows: 



. . e x x k f e\ k 1 
*fcW < nr-r T = 7k ex P 



. py^- (x + pv^+ijlogfl + ^z 
exp[£(p,x)] (B.5) 



V 2 



7TX 



Now, keeping x fixed, we consider three different regimes for k (or equivalently for p). 

First, for x < k < ^ (equivalently < p < we can expand the logarithm as an absolutely 
convergent series, and collect powers of x" 1 / 2 : 

E(p,x) = P^-{ X +P^+l){^-^ + ^j- 2 -^ + . 

P 2 , a t(p) , n ^ 

= -y + L^ 

i=i 
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where the coefficients ai{p) are polynomials in p: 

Oi(p) = R i+ V 

Hence 



p 2 1 



(« + !)(» + 2) 2i 



(B.7) 



where the &i(p) are also polynomials in p which we do not give explicitly. 

Second, for ^ < k < e 2 x (equivalently < p < (e 2 — l) v / x), we use the fact that tfe(x) is 
monotonic decreasing in k (for k > x and x fixed), to give 

**(*) < t ¥ {x) = -^y|exp|Q-|log|)x} = t (2) (x) (B.9) 

which is (clearly) p-independent. The coefficient of x in the exponent is negative. 
Finally, for k > e 2 x (equivalently p > (e 2 — l)y / x), we have 

,og(l + -|)>2 (B.10) 

and hence 

t k (x) < -—e"^ = t {3) (k,x) (B.11) 
\/2irx 

in which we emphasise the p-dependence, since we will be integrating over p. 
B.3 A lower bound for g^ m (x) 

Using the previous definitions of g n ( x ) and tk(x) (eqs. (|4.3p and (|B.3[) respectively), we have (for 
n > x) 

n oo 

<7«»(s) = £t fe (x) = 1- E 

fc=0 k=n+l 

1 - /,f /2 (k, x)dk - f£j 2 1 {2) (x)dk - f°? x t {3) {k, x)dk {x<n< 3x/2) 

1 " fn 2x t {2 \x)dk - f^ x t {3) (k,x)dk (3x/2 < n < e 2 x) (B.12) 

l-/~f (3) (fc,:c)dfc (e 2 x<n) 



> 



coo re x poo 

> 1-1 t {1 \k,x)dk- / t (2 \x)dk- / t {3) (k,x)dk 



3x/2 



(B.13) 



where the t^\k,x) are the upper bounds for tk(x) derived in Appendix IB.2I We have also used the 
fact tk(x) is monotonic decreasing as a function of k. 

In the same way that we introduced p as a scaled proxy for k (eq. (|B.4|) ). we similarly introduce 
m as a scaled proxy for n: 

n — x ,„ 
m = — =- B.14 

\/x 
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The first integral can then be written 



OO / \ 

= l"*M + £^ (B.15) 

i=l 

where $(m) is the cumulative normal function, and the Cj(m) have no x-dependence. 

The second and third integrals have the forms ax 1 / 2 e~@ x and r yx~ 1 ^ 2 e~ Sx respectively, where a, 
(3, 7 and 5 are constants: 



a = 


e 2 - 2 

e 2 




V37T 




P = 


^ l0g ( 


1) 


7 = 


1 








5 = 


1 + e 2 


> 



(B.16) 

The important thing to note is that these two integrals are exponentially small when compared with 
the first integral. 

Hence, we can write (as an exact result) 

OO / \ — fay. 

i=l * 

This is true for any m > (equivalently, for any n > a;). However, we see that the right-hand side has 
a well-defined asymptotic limit as x — > oo only for fixed m (i.e. and not for fixed n). 
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